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Abstract. An explicit perturbative computation of the Mori’s memory function was performed by Gotze 
and Wolfle (GW) to calculate Generalized Drude scattering (GDS) rate for the case of electron-impurity 
and electron-phonon scattering in metals by assuming constant electronic density of states at the Fermi 
energy. In the present investigation, we go beyond this assumption and extend the GW formalism to the 
case in which there is a gap around the Fermi surface in electron density of states. The resulting GDS is 
compared with a recent one by Sharapov and Garbotte (SC) obtained through a different route. We find 
good agreement between the two at finite frequencies. However, we hnd discrepancies in the dc scattering 
rate. These are due to a crucial assumption made in SC namely uj » \S{€ + lj) — X'*(£)|. No such high 
frequency assumption is made in the memory function based technique. 

PACS. 74.25.Gz Optical properties - 72.10.-d Theory of electronic transport; scattering mechanisms 

1 Introduction 

The study of transport properties like optical conductiv¬ 
ity is very important to understand the electronic inter¬ 
actions in complex many body systems like cuprates pQH]. 

The electronic interactions comprises of electron-phonon, 
electron-boson (spin-fluctuations), electron-impurity, elec¬ 
tron-electron interactions. Experimentally, the signatures 
of these interactions can be grasped by using optical data 
((t(w, r))I31ll] which includes the deduction of the Gen¬ 
eralized Drude scattering (GDS) rate and mass enhance¬ 
ment factor using the standard form 

UJp 1 

^ ^ Att 1/t{uj,T) + iuj^l + X{uj,T)) ^ ^ 

Here 1/t{uj,T) is the frequency and temperature depen¬ 
dent scattering rate, A(a;, T) is the frequency and temper¬ 
ature dependent mass enhancement factor and ujp is the 
plasma frequency. On theory side, the derivation of the an¬ 
alytical formulae for these quantities (like 1/t{uj,T) and 
X{u!,T)) is very complicated. 

Generally, the optical conductivity is calculated us¬ 
ing Boltzmann’s equation by assuming relaxation time 
approximation[5l[6l[7]. In that picture, scattering rate is 
considered as constant (independent of frequency and tem¬ 
perature) and the resulting behaviour corresponds to the 
Drude behaviour. In systems such as cuprates where the 
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electron-boson, electron-electron interactions are impor¬ 
tant, this approach is inadequate as experiments show that 
scattering rate does depend on frequency and temperature 
[8l|9]. The complete solution of the transport problem in 
cuprates and other strongly correlated materials is compli¬ 
cated as there are no controlled perturbation parameters. 
For example, in simple metals where, as first shown by 
Holstein[TU], the perturbation parameter is Vg/vp (sound 
speed/Fermi velocity) which is a small parameter and per¬ 
turbative calculations are justified. Building upon Hol¬ 
stein’s work on metals[Tn], Allen[TT] has derived a relation 
for frequency dependent scattering rate. This has been 
further generalized by Mitrovic and Fiorucci[12j by con¬ 
sidering the effects of non-constant density of states. Fur¬ 
ther this has been extended recently for the finite tem¬ 
perature case by Sharapov and Garbotte[T3]. All these 
approaches are based on some assumptions such as ne¬ 
glecting vertex corrections. To go beyond this assump¬ 
tion, we have developed a formula for frequency and tem¬ 
perature dependent scattering rate using memory func¬ 
tion techniaue |14l[T5] which includes the effect of vertex 
corrections [in] (where the current-current correlators are 
directly computed without writing them in terms of single¬ 
particle Green’s function). This technique is a general¬ 
ization of Zwanzig projection operator techniaue |171ll8) . 
Physically, this approach is very appealing, because the 
conductivity < 7 ( 0 ;, T) can be cast into the generalized Drude 
form with frequency and temperature dependent scatter¬ 
ing rate. Recently, it has also been used by several authors 
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to study the transport properties of different systems [I^ 

The paper is organized as follows. In Sec.II, we have 
elaborated the Gotze-Wdlfle memory function formalism 
m- In Sec.Ill, we go beyond the constant electronic den¬ 
sity of states assumption and introduce the gapped density 
of states and calculate the imaginary part of memory func¬ 
tion. Here we also discuss the dc and ac imaginary part 
of memory function in different temperature regimes and 
in appropriate limits we reproduce GW results. In Sec.IV, 
we compare our findings with SG result [15] and finally we 
conclude with a brief discussion in Sec.V. 


2 Gotze-Wolfle Formalism for 
electron-phonon scattering 

In this section, a short introduction to GW formalism is 
presented[TB]. The Hamiltonian used for electron-phonon 
interaction is given by 

H = Ho + H,p + Hp, ( 2 . 1 ) 

where Hq is the Hamiltonian for non-interacting electrons 
and is represented as 

Hq = c(l^)Ck^cr''k,(T) (2-2) 

kjtr 

where e(k) is the band dispersion and c|^^, Ck.o- are cre¬ 
ation and annihilation operators with wave vector k and 
spin a. The Hamiltonian Llep represents the electron - 
phonon interaction and is given by 


where J = en(k)c|^ the current density and 

k 

z;(k) is the velocity dispersion. Here [J{t),J] denotes the 
commutator, (...) denotes the ensemble average at tem¬ 
perature T and ((...)) denotes the Laplace transform of 
the ensemble average. 

According to the Gotze and Wolfle approach[T6j, the 
memory function is defined as 


M{z) = z 


X{z) 


Xo - x{z) 


or 


Xiz) = Xo- 


M{z) 


■M{zy 

(2.7) 

where xo corresponds to the static limit of correlation 
function (i.e. x o = He/m) [ig. Using this, the conduc¬ 
tivity from eqn.(2.5l in terms of memory function can be 
written as 


r(z) = 




Att z M{z) 


( 2 . 8 ) 


In ref.[T6|, an expansion for M(z) = 1 1 + 

Xo V Xo 

is used. Basis of this assumption is the smallness of electron- 
phonon interaction energy as compared to the Fermi en¬ 
ergy of free electrons. Using this expansion and on keeping 
the leading order term, the memory function M(z) can be 
written as 


M{z) = z 


xiz) {{J;J)). 


= z- 


(2.9) 


Xo Xo 

To compute memory function, we need ((J; J))z which by 
using equation of motion is 


z((J; J))z = ([J, J]) + (([J, Hep]; J))z- (2.10) 


Hep = ^ |^Z)(k - k')C|{^, ^Ck.cT&k-k'+ h.c. 

k.k' .(T 


(2.3) 


Here &k-k' j ^k-k' annihilation and creation opera¬ 

tors for phonons and H(k —k^) is the electron-phonon ma¬ 
trix element. Here the symbol h.c. corresponds to the Her- 
mitian conjugate of first term. The third part of eqn.(2.11 
represents the free phonon Hamiltonian 




1 + - 


(2.4) 


where uiq is the phonon frequency. 

According to the linear response theory, the dynamical 
conductivity is defined as[i 




(2.5) 


Here u/^ = AirNee^/m is the square of plasma frequency 
where e electronic charge, m electron mass and Ng is the 
electron density, z is the complex frequency and x(^) is 
the current-current correlation function defined as 

/*00 

x{z) = {{J-,J))z=i e^^\[J{t),J]), (2.6) 

Jo 


As the first term of r.h.s is zero, hence the above expres¬ 
sion is equivalent to second term which can be further 
calculated by applying equation of motion. 


z{{[J,Hep];J))z = ([[J,Hep], J]) - (([J, Hep]; [J, Hep])).. 

( 2 . 11 ) 

For z = 0, ([[J,Hep], J]) = (([J,Hep]; [J, Hep])).=o. Thus, 
the memory function M{z) becomes 


M{z) = 


(/)(0) - (/{z) 

zxo 


( 2 . 12 ) 


Here (j){z) (called as correlation function) is defined as 


,^(z) = (([J,Hep];[J,Hep]))^. (2.13) 


For the present case of electron-phonon interaction, the 
correlation function from eqns.(2.3) and (|2.13|) is 


bi(p) - ^'l(p')] 

k,k' p,p' O',dr' 

X ((H(k - k')c^Ck'.rT^k-k' - D*{k- ^)b\^_^,c\^, 

D{v - P')Cpy Cp,_^r6p_pr - H*(p - P')^p-P'Cp'.T w)) 

(2.14) 
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XI K(k) - (k')] bi(p) - 1^1 (P')] 

k,k' p,p' cr.cr' 

X [i^(k-k')i?*(p-p')] 

X (((cXck'..&k^k';^>l,-pTj,,^^,Cpy)) 

+ ((ck,.4',aCk';^-p'S'.<r'cX'))) ■ 

(2.15) 


To evaluate the (('(z), we need to calculate 
,. kt J 

P-P' p ' a P.tr 

lated as (using definition 2.6) 


(('^k o-'^khtT^k-kd _ /c', iC„„i)) which Can be calcu- 


((4,<TCk'.,T&k-k'; &p_pTp,_^,Cp = 

[/(I - /')(1 + n) - f (1 - f)n] 


2 - ek' + ek - Wk_k' 


Here / = /(ck) = + l) ^ and n = n{uju-k') = 

(g^“k-k' _ represent the Fermi and Bose distribution 
functions and /? corresponds to inv erse of temperatu re. In- 
serting this equation in eqn.( 2.151 and hence in eqn.( 2.121 


M'\oj,T) = 


2tt 1 
3 mNf 


^|l?(k-k')P(k-k')Y(l-/)r 


k,k' 


- 1 


to 


<5(ek-ek' -Wk_k' +w) 


+ (terms with a; —>■ —w)]. 


2 r°° r°° 


/•oo ptt 

/ dk'k'^ / (i0sin0(5((7 — |k — k^) 
do do 


/'(I - f)n 


- 1 


<^(ek - Ck' - tuk_k' + w) 


+(terms with a; —>■ —w)]. 


Here due to the presence of Fermi factors /'(I — /) the 
integrand has finite value only around the Fermi surface 
and vanishes outside the strip of width 2//3 (cu << ep)- 
Thus k and k^ can be approximately replaced by kj?. With 

r 

this the 9 integral can be simplified as / d9 sin 96 {q — 

do 

V2kFVl — cos 9) which will yield the result y^. Using 

kp 

this and converting k integrals into energy integrals, the 
above equation reduces to 


M''{uj,T) = 


4 N^m^ep 
3 {2pfN,kl Jo 


CQD 


dqq^\D{q)\^^ 


((4.aCk'.a^k-k';&p_p'C^, ,^'Cp ,^/ 

POO 

i / dte***([cj,^(t)ck',a(0^k-k'(0 

;^p-p'i'ys,-'])- (2.16) 

Using Ck,(T(t) = Ck,o-e“*'^‘‘‘ and performing the integration 
over time and ensemble average, we have 


de 


n 


-oo 

1 - 1 

g/3(e-ej7+a;-a;,) _|_ y 

+(terms withw —>■ —w)]. 


( 2 . 20 ) 


This is an expression for the imaginary part of memory 
function as deduced by Gbtze-Wolfie[l6] . It can be simpli¬ 
fied by using electron-phonon matrix element for acoustic 
phonons which is defined as[S] 


.(2.17) 


7d(q) = 


1 


2miNojo 


- 1/2 


qC{q); w, = c^g, (2.21) 


where C{q) is the slowly varying function of q, rrii is the 
ion mass, N is the total number of unit cells and c, is the 


sound velocity. To analyze the eqn.(2.20), various limiting 


and then by taking the limit z —>■ w -|- i?/, ly —?> 0^, the 
imaginary part of the memory function can be expressed 
as 


cases using eqn.(2.21) were discussed in ref. [16] . 


(2.18) 


Convert the summations over k and k^ into integrations 
and assuming that k is pointing along the z-direction and 
k subtends an angle 9 with it (at the end k integration 
over all directions and magnitudes is to be performed). 

Insert an integral dqd{q— |k — k^) over q to stratify the 
calculation as given below. Thus the eqn.(2.18l becomes 


3 Memory Function with Gapped Density of 
states 


In this section we go beyond the assumption of constant 
electronic density of states and we consider a system with 
a gap around the Fermi surface. In this case, density of 
states is zero in energy region {—A, A). Thus the energy 
integration in eqn.(2.20) has to be 


p^fp — 


I = 


' —oo 
/•oo 


de 


de- 


pl3{e-eF) 


1 


g^(6-ej^) I ^0{e-eF-\-(A-Uq) _j_ ^ 

g/3(e-ej^) ^ 


Jep+A + I g/3(e-eF-Hc^-i^,) + 1 ' 

After simplification we have 

7=1 ' 


(3.1) 


(2.19) 


^ g/3(i^-w,) _ I 

/ (1-7 e^(-‘^+“-‘^‘>))(l-7 e^^)e^(‘^- 
V (1 g/3(zi-ha;-a;,))(l g-/3zi) 


(3.2) 
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Using this, the imaginary part of memory function can be 
written as 


M''{uj,T) = 


4mkp 
- 1 


/o 


I 

dqq^\D{q)\'^-n 

1 


ui - 1 

1 _(_ e/3-4 \ / ^ ^-P(A-uj+ujg) 


p_ /■‘^D 

M"{oj,T) = — dna^F{n) 
^ Jo 


- 1 


1 


log 


g/3(w-t2) _ 
2 _|_ g—/3(^—<^+^2) ' 


g/3t2 — I 2 _j_ g— 0(A-\-uj —12) 

— (terms with a; —>■ —w)], 


where a^F{n) is defined as 


a‘^F{f2) = 






(3.4) 


(3.5) 




This is known as phonon spectral function|ll]. In the case 
of cuprates, it is replaced by /^x(f?) which represents the 
boson spectral function|13j. This form is same as given by 

Allen m (a^Fin) = ^((| (t/(k) - t;(k'))" 

6{hf2Q - hn)))). 

Equation (3.4) is our main result. To discuss it in 


various temperature and fr equen cy regimes, we use the 
phonon matrix element eqn.(2.21) and calculate M"{uj, T) 
in next subsections. 


3.1 DC memory function 

In the zero frequency limit and assuming C{q) as a con¬ 
stant i.e. C{q) = I/pf [S], the imaginary part of the mem¬ 
ory function (eqn.(3.3)) becomes 


1 N 

o riirriiUp Jq 


— log 


Wo 


^g/3co, _ ]^^^g /3oj, _ 

I-be^^ 1 _|_ g-/3(-4-|-:^,)- 


I -|- 


(3.6) 


Now consider the case of T » ujd,^, the above equation 
reduces to 


M"(0,r) = -r 


N 


rto 1 -1. 

dqq^ 


■ f3Wg 


^ ' 1 + g-/3zi y g/3-4 

-|-(terms withw —> —w)]. (3.3) 

This is the desired expression for the frequency and tem¬ 
perature dependent imaginary part of memory function. 
For A = 0 and using phonon matrix element (eqn.(2.21)), 
this expression reduces to the expression (refer eqn.(54(a))) 
given in original GW work[T6], as it should. In actual prac¬ 
tise (i.e. for an arbitrary form of gap around the Fermi sur¬ 
face) , the general expression of the imaginary part of mem¬ 
ory function is complicated and is difficult to proceed ana¬ 
lytically. A general formulae is given in appendix ([A|). Thus 
for the simplicity of calculation, we have discussed the spe¬ 
cific system in this article. Further to write M''{uj,T) in 
compact form, change the variable coq to 17 in above equa¬ 
tion which can be rewritten as 


X log 

On substituting x = 


8 mmik% Jq tUq (/Swq)^ 

2-j3A-j3ujq' 


2 — /3A -|- Jujq _ 
qOo 


(3.7) 


(i.e. Pojq = x) where ©f is the 


qoT 

Debye temperature, the dc memory function reduces to 


M"(0,r) = -TT 


1^3 N 6 
8 ^ mmik^OD^^yODj 


np^D 

/ dxx'^ix — l)log 
Jo 


2- pA-x 
2- pA + x 


(3.8) 


This expression under case T » ujj^, A is equivalent to 


M"(0,T):xA(^ +A+ 1 


A2 


8A Od 

Od ' Od ' T \Wd ' ^ fT 


(3.9) 

where A refers for constant numerical factor. Similarly for 
T « u}d,A, the eqn.(3.6) becomes 


^ s'" mm,k%OD^^{o 


T Y 


rP0n 


dxx 


.4„-a 


log 




This expression can also be simplified as 


M"(0,T) ~ 


1 


5 4 XOd 


T 


.(3.10) 


(3.11) 


Substituting the equation (3.10) in eqn. (2.7) and hence in 


eqn.(2.5), leads to the expression of dc conductivity for the 


electron-phonon interaction. Here if we insert gap A = 0 
in the eqn.(3.6), we obtain eqn.(54(b)) as given in ref. [TH] . 
as expected. 


3.2 AC memory function 


We proceed again with eqn. (3.3) to study frequency de¬ 
pendent behaviour of memory function in different regimes. 
In the high frequency regime i.e. for w >> ujd and using 
same approximation (C{q) = 1/pf) as considered for the 
dc case, the imaginary part of memory function becomes 


I 


M''{UJ,T) = -TT' 


N 


^ 5 I n 

dm 1^ - 


8 mmikp Jq Pujq uj 


log 


I -b e 


-PiA-u,) 


1 -|- e“/3(-4+‘-.') 




log 


I -b 

I -I- q-P{A-u>) 


(3.12) 
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When the gap is smaller than the \oj — ujd\ i-e. Z\ < |a; — 
wul, the above equation reduces to 0 


M'\uj,T) = -r 


N 


-<lD 


T 


D 


8 mmikp0D 

J dxx"^ coih y—j . (3.13) 

From this we identify that at high temperature, the imag¬ 
inary part of memory function becomes temperature and 
frequency independent. This means the saturation behavio¬ 
ur of M"{(jj,T) for w >> ujd- The reason is that under 
this condition, the integral approaches to (0n /T)^ and 
it cancels with prefactor {T/0p,)^ in eqn.(3.13). At low 


temperature, it varies linearly with temperature as the 
integral approaches to {0d/T)‘^ 


In the next section we compare our findings (eqn.(3.4)) 
with the Sharapov-Carbotte|13j. 





<a(eV) 

(b) 



4 Comparison with Sharapov-Carbotte results 

Sharapov and Carbotte has deduced a relation for the gen¬ 
eralized Drude scattering rate [13] taking electron-boson 
interaction and non constant electronic density of states. 
Using Kubo formula|32] and calculating the self energy 
under certain assumptions (as discussed below), they de¬ 
rived the following expression 


1 


"(w,r) 


- / dQI^xi^) / duj' 

^ Jo J — oo 

N{uj'-n) ^ N{-uj' + n) 

iv(o) iv(o) 

X [n(l7) -I- f{Q — w')] 

X [/(w' - w) - f{uj' +u })], 


(4.1) 


In the opposite case — ujd\ < A, eqn.(3.3l leads to van¬ 
ishing scattering rate. 


Fig. 4.1. Comparison plot of scattering rate 
1/t{ijj,T){=M"{ijj,T)) calculated using Memory function 
approach (Mem, solid, green) and Sharapov-Carbotte ap¬ 
proach (SC, dashed, black) at temperature T= lOK and 
200K and at gap A = 0.02eV and 0.20eV. The agreement is 
excellent. 


0.010 

0.008 


0.004 

0.002 


W'0,A=0.02 / 

0.025 

0.020 

> 0.015 

w=0.05,A=0.02 / 




< 0.010 



0.005 

0.000 



20 40 60 80 100 


20 40 60 80 100 


T(K) 

(a) 


T(K) 

(b) 


where /^x(^) is the boson spectral function and N{uj) is 
the quasiparticle electronic density of states and iV(0) is 
for normalization. In deriving the above formula, the fol¬ 
lowing assumptions were made: ( 1 ) vertex corrections were 
neglected, ( 2 ) energy independent character of plasma fre¬ 
quency in the vicinity of Fermi level and (3) \E{e + oj) — 
A'*(e)| « OJ where S{e) is the ele ctro nic self energy. 

T o com pare our approach (eqn.( 3.4)) with SC approach 
(eqn.( |4.l[ )), we have done calculations using models for 
electronic density of states and the boson spectral func¬ 
tion. First in SC approach, for the electronic density of 
states, we use a square well type model with center at 
Fermi energy and considered a gap of 2Z\ around it. Same 
gap is taken in our approach (eqn.(3.4)). Second, for the 
boson spectral function, /^x(l7), we modelled it as Lorentz- 

rn 

ian of the type 777 - . „ - , , „ where f2E represent the 

(sj - iiEr -h (I r 

boson peak frequency and F is the width of the Lorentzian 
(this form has been used extensively in ref. [551IMII55] 1. 



0 20 40 60 80 100 

T(K) 

(c) 

Fig. 4.2. Temperature variation of scattering rate with two 
different approaches namely memory function (Mem, dotted) 
and Sharapov-Carbotte (SC, dashed) at gap 0.02eV. (a) dc 
case (b) at a; = 0.05eV (c) at a; = O.SeV. 


Thus for comparison, we use the same form of /^x(i7) 
in SC approach and our approach. In the whole analy¬ 
sis, we have fixed the value of Qe and F as 0.02eV and 
0.04eV respectively in both approaches. The value of De¬ 
bye frequency (the upper limit of phonon frequency in¬ 
tegration eqn.(3.4)) is very much high as compared to 
the Lorentzian width, hence ojd doesnot give any effect 
in whole calculation. To compare the results from both 
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the approaches, the frequency dependent scattering rat e 
has been plotted at different temperatures. In fig.(4.1), 
we can observe an excellent agreement between both the 
approaches. As the gap magnitude is increased, the scat¬ 
tering rate shows suppressio n upto the frequen cy a; ^ Z\ 
as expected (compare figures [T 1(a) and 4.1(c)). These re¬ 
sults are qualitatively agrees with experimental results [HI 

m- 


In fig.(4.21, we plot l/r(a; —)■ 0,T) as a function of 
temperature T. Here we can observe that the memory 
function approac h yield s more magnitude over the SC 
approach. In fig.( |4.2(^ ) i.e. in zero frequency limit, the 

I/tmf - 17^^ ^ andl/rsc rep¬ 


ratio 




lOOK 


resents the scattering rate by memory function technique 
and SC technique respectively, is 0.7 which becomes 0.4 
at w = O.OSeV (as shown in fig.(4.2(b))) and at a; = 0.5eV 
it further reduce to 0.031 (as shown in fig.(4.2(c)). This 
shows that the difference between scattering rates using 
memory function approach and SC approach reduces as 
we go from dc limit to finite frequency limit. Also both 
approaches explain the Holstein ’ s mech anism at T = OK 
[^[37] (as shown in fig. ( |4.2(b) 4.2(c))). Thus we notice 
that there are discrepancies between the two approaches 
in the d.c. limit. The reasons are discussed in next section. 



increase with increasing temperature. These discrepencies 
observed in the dc limit are discussed in the next section. 


5 Discussion 


It is observed that the finite frequency scattering rate us¬ 
ing memory function formula is in excellent agreement 
with the same obtained from SC formula as shown in 


fig. (4.1). This shows that the assumptions made in these 


two different approaches are consistent at finite frequen¬ 
cies. However, while discussing dc scattering rate, we ob¬ 
serve significant discrepancy between the two approaches 
(fig.(4.2(a) I and fig.(4.3)). To illustrate it further, we have 
plotted the difference in the magnitudes of scattering rates 
calculated by both approaches. The di fferen ce (1 /tmp’— 
1/t'Sc) at a = 0.02eV is plotted in fig.(5.1). Here we find 
that this difference increases with the rise of temperature. 
The reason behind this difference in the dc case is the as¬ 
sumption made by SC i.e. w >> lA'(e-l-CLi) — A'*(e)| which 
becomes more severe in high temperature regime. To clar¬ 
ify this fact, we plot the quant ity \ S{e -I- w) — A'*(e)| as a 
function of temperature in fig.(5.2) (where the expression 
used for S{uj) has been given in ret. [TH] 1. It shows that the 
magnitude of the difference of self energy increases with 
the temperature. This shows the stronger violation of the 
condition w » 127(6 -I- w) — A'*(e)| in high temperature 
limit. It implies that SC formalism is not appropriate to 
study the dc behavior and the disagreement is severe at 
high temperature, but it is quite reasonable for the finite 
frequency case. 

Regime of validity of Memory function: It is important 


0.008 

0.006 

0.004 

0.002 

0.000 


A=0.02 


0 20 40 60 80 

T(K) 


100 


Fig. 5.1. Variation of difference = (1/tmf — I/fsc) of dc 
scattering rate with temperature calculated by two different 
approaches MF and SC at A = 0.02eV. 


Fig. 4.3. Comparison of dc scattering rate (l/r^c) as a func¬ 
tion of A using Memory function approach (Mem, Purple) and 
Sharapov-Carbotte approach (SC, Red) at various tempera¬ 
tures 50K, lOOK and 200K. 


Next, we have plotted l/r^c at different temperatures 
as a function of A and compare the both approaches. Here 
we observe that decreases with the increase of gap 

energy A. Also, we find that the difference between the 
l/rsc and 1 /tmf is not much dependent on A, but it does 


to recognize that memory function calculation of l/r(a;) 
presented here is also perturbative in character and has 
its own limitations. Here in calculation, we have used the 
expansion for M{z) which is given below eqn 2.8 where 
M{z) = {zxiz)/xo) (1 + x{z)/xo + ••■)• In the leading or¬ 
der approximation, we took xiz)/xo << Ij which im¬ 
plies that the magnitude of memory function is smaller 
than z {M{z)/z = x{z)/xo)- Thus to ensure this point 
we have plotted the frequency dependent scattering rate 
for temperature 10 K and 200 K at A = 0.02eV and 
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Fig. 5.2. Plot of \S{e + ca) — L'*(e)| with temperature at dif¬ 
ferent frequencies such as (a) uj = O.OOOleV, (b) uj = O.OOleV 
and (c) u) = O.OleV . Here S{u)) represents the self energy and 
* corresponds to the conjugate. 


A Appendix-1 

The imaginary part of memory function (eqn |2.1^ is 
M"(ca,r) = \Dik-k')\^k-krf{l - f)n 

® k,k' 

■ e‘a/T _ 2^ 

-(5(e -e' - Wk_k' + w) 

UJ 

-|-(terms withca —>• —w)]. (A.l) 

Converting the summations into energy integrals and in¬ 
serting the dq integral as before, we have 





>ii(eV) <a(eV) 


(c) 


(d) 


Fig. 5.3. Solid line: Variation of scattering rate with frequency 
at different temperature and A (in eV). Dotted line: Plot of 
f(uj) = UJ to check the validity of memory function approach. 


M''{uj,T) = 


27r 

3 mNe 7q 


/»CXJ ^OO /»CXJ 

/ dq deN{e) / ddN{e') 

Jo J — oo J — oo 

[ dOsineSiq- \k - k'\)\D{q)\'^q'^ f {1 - f)n 
Jo 

_ I 


UJ 


-5{e — e' — ujq + uj) 


-|-(terms with a; —>■ —w)]. 


(A.2) 


Here the energy dependent density of states N{e) has been 
introduced. Thus on solving the integrals over e' and 0, the 
above equation reduces to 


M"{uj,T) = 


2 tt iV 2 r 

3 mNekpUj Jq 


deN{e) 


N(^e — ujq uj) 


dq\D{q)W^ 

g^(e-eF) I 

ed(e-er-) -(- 1 el^‘^9 — 1 

- 1 


g/3{e-ep-uiq+ui) _|_ ^ 


— N{e — UJq — Uj) 


gP(€->LF-UJq-U>) _|_ 1 


(A.3) 


0.20eV in fig.(5.31 an d com pared with the linear variation 
f(uj) = UJ. From fig.(5.3), we find that our app roach i s 
valid for regime where A » T (as shown in fig.(5.3(a)), 


5.3(c) and 5.3(d)). Thus it is quite suited to study the 
pseudogap phase of cuprates when A is greater than the T. 
But our approach is not good in the low frequency regime 
to discuss the case where Z\ ~ T (as shown in fig.( 5.3(b)| )). 
Such small gap scenarios occur in spin/charge density 


This is the general expression for the imaginary part of 
memory function (called as CDS). 
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